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ABSTRACT 

Context. The accretion histories of embedded protostars are an integral part of descriptions of their physical and chemical evolution. 
In particular, are the accretion rates smoothly declining from the earlier toward later stages or in fact characterized by variations such 
as intermittent bursts? 

Aims. We aim to characterize the impact of possible accretion variations in a sample of embedded protostars by measuring the size 
of the inner regions of their envelopes where CO is sublimated and relate those to their temperature profiles dictated by their current 
luminosities. 

Methods. Using observations from the Submillimeter Array we measure the extents of the emission from the C**0 isotopologue 
toward 16 deeply embedded protostars. We compare these measurements to the predicted extent of the emission given the current 
luminosities of the sources through dust and line radiative transfer calculations. 

Results. Eight out of sixteen sources show more extended C**0 emission than predicted by the models. The modeling shows that the 
likely culprit for these signatures is sublimation due to increases in luminosities of the sources by about a factor five or more during 
the recent 10,000 years - the time it takes for CO to freeze-out again on dust grains. Eor four of those sources the increase would 
have had to have been a factor 10 or more. The compact emission seen toward the other half of the sample suggests that C'*0 only 
sublimates when the temperature exceeds 30 K - as one would expect if CO is mixed with H 2 O in the grain ice-mantles. 
Conclusions. The small-number statistics from this survey suggest that protostars undergo significant bursts about once every 
20,000 years. This also illustrates the importance of taking the physical evolutionary histories into account for descriptions of the 
chemical structures of embedded protostars. 
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1. Introduction 


In their earliest stages, solar-type protostars are characterized by 
large amounts of cold gas and dust surrounding them. As the 
young star begins to heat up and accretion proceeds, this en¬ 
velope is dissipated revealing the pre-main sequence star sur¬ 
rounded by a circumstellar disk. Studies of young stars during 
the embedded stages are of prime importance for our understand¬ 
ing of their formation and early physical and chemical evolution: 
this is the phase when the star accretes the bulk of its mass and 
a circumstellar disk is formed. The physical and chemical struc¬ 
ture of this disk possibly sets the initial conditions for subsequent 
planet formation. One of the important questions for our theories 
of formation of stars is how accretion proceeds: is it taking place 
at a relatively constant rate throughout the protosteiiar evolution 
or strongly varying, e.g., related to the formation of the disk and 
possible fragmentation within it? Indirect evidence suggests that 
the accretion rates may be varying through the protostars’ early 
years - for example observations of outflows that in some cases 
show clear knots or bullets that can be attributed to changes in 
the un derlying accretion rates (e.g., Reipurth||1989[ Arce et al. 
2013 I. This paper presents maps of the C**0 emission on scales 
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of a few hundred to 1000 AU toward the dense envelopes of a 
sample of 16 embedded protostars with the Submillimeter Ar¬ 
ray (SMA). The aim of this work is to address whether possible 
variations in the protosteiiar accretion rates can be traced by their 
chemical signatures. 


Strong luminosity bursts have for some time been recognized 
as one of the characteristics of young stars - at least during spe¬ 
cific times of their pre-main sequence phase where they may ap¬ 
pear as FU Orionis objects at visible wavelengths (see, e.g., |Au-| 
|dard et al.|2014[ for a recent review). Such outbursts can be tied 
to the accretion onto and through the central protosteiiar disks 


- and possible instabilities in those (e.g.,|Bell & Lm|l994| Ar-| 

mitage et al. 

20011 Vorobyov & Basu||20051 |Zhu et al. 2009 

Martin & Lu 

5Ow||2011|). Due to the significant amounts of ob- 


scuring gas and dust, it is still unclear what happens in the ear¬ 
lier protosteiiar stages. Photometric comparisons of nearby star¬ 
forming regions at different epochs show small-scale luminosity 
variations on time-scales ranging from days to years (e.g., |RebuII] 
et al.|2014|) - and are able to identify individual sources undergo¬ 


ing bursts in intervals of 5000-50,000 yr ( |Scholz et al.|2013T ). If 
these bursts are an indication that the accretion rates are indeed 
also strongly varying, it would definitely be a critical aspect of 
the protosteiiar evolution. 
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In the earliest stages the luminosity of the protostar is 
thought to be dominated by the release of gravitational energy 
due to the accretion onto the central star (L ~ GM^MIr^). Gener¬ 
ally, protostars are found to show a very broad luminosity distri¬ 
bution and be under-luminous compared to the expectation from 
constant infall toward the center of their natal cores (see, e.g., 
Kenyon et al.|| 199^ and [Evans et al.||2009[ as well as [Dunham] 


et al.|2014 for a recent review). This is another piece of evidence 


that has been taken in favor of episodic accretion during the ear¬ 
liest protostellar phases. In that picture, protostars spend most of 
their lifetimes in a low accretion rate mode, assembling material 
in the disk, and thus appear less luminous than what should be 
expected based on accretion from their envelopes directly onto 
the central star. They could then accrete significant mass in rel¬ 
atively short bursts accompanied by increases in their luminosi¬ 
ties that statistically would be relatively difficult to observe. Still, 
because of the complex environments of young stars and possi¬ 
ble continued infall onto and through the disk also during the 
more evolved stages, this picture may be too simplistic and the 
luminosity distribution itself not offering a good constraint (e.g. 
[Padoan et al.|201^ . 

Besides the obvious important dynamical consequences of 
the exact protostellar accretion histories, they may also have im¬ 
portant consequences on other aspects of the early evolution of 
protostars - for example, the chemical processes leading to the 
formation of larger molecules. This process involves a complex 
balance between molecules freezing out onto dust grains, grain- 
surface chemistry and eventual evaporation of the species before 
they may be (re-)incotporated into ices in the circumstellar disks 
(see, e.g., [Herbst & van Dishoeck|[2()09[ [Caselli & Ceccarelli| 
|2012| for recent reviews). Variations in the accretion rates and 
consequently luminosities of the protostars will affect how long 
given molecules are present in solid or gaseous form and what 
the time-scales are for the regulating chemical processes at dif¬ 
ferent points during the protostellar evolution. 

An interesting question is how these two issues, the accre¬ 
tion histories of protostars and corresponding luminosity evolu¬ 
tion on the one hand and chemical evolution on the other, cor¬ 
respond. The strong dependence of the protostellar luminosity 
on the accretion rate will reflect directly in the thermal struc¬ 
ture of the envelopes through the (almost instantaneous) heat¬ 
ing of the dust ( Johnstone et al.|2013[ ). As these times are also 
characterized by freeze-out and sublimation of molecules that 
are strongly temperature dependent, clues may be found in the 
composition of ices as well as distributions of molecules in the 
gas- and solid -phase (e.g.,|Lee|2007 Poteet et al.|2013 Visser & 


|Bergin|2012[ [Visser et al.|2015[ ). One example of this is obser¬ 

vations that show significant amounts of pure CO 2 ices present 
around low to moderate luminosity protostars ( Kim et al.|2012| l. 
The presence of pure CO 2 ice is a good indication that significant 
thermal processing of the ices has taken place prior to the cur¬ 
rent evolutionary stage of the protostars. Still, such observations 
of the total column density of a specific ice species only reveal 
the integrated history of the source - and cannot be used as con¬ 
straints on whether this distillation has taken place during the 
protostellar evolution or perhaps prior to the onset of collapse. 

Another way of addressing this issue may be through re¬ 
solved observations of the distributions of molecules in the envi¬ 
ronments. One possible example of this are ALMA observations 
of the deeply embedded protostar IRAS 15398-3359 pprgensenj 
jet al.|20l3|) . The resolved images show a depression in the emis- 
sion of H‘^CO^ in the inner 150-200 AU of the central proto¬ 
stellar core. One explanation for this depression is that HCO^ is 
destroyed by reactions with extended water vapor. However, the 
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HCO^ depression is seen on the scales where the temperature 
is as low as 30 K given the current luminosity of the protostar 
- well below the temperature of 90-100 K needed for water to 
sublimate. A possibility is that the source has undergone a burst 
in accretion during the last 100-1000 years, increasing the lumi¬ 
nosity by up to a factor 100 above what it currently is, causing 
water to sublimate on much larger scales. Ideally, one would use 
similar observations for many sources to constrain the promi¬ 
nence of such possible bursts: if such variations are indeed due 
to accretion variations, the distributions of even simpler species 
such as the optically thin isotopologues of CO may in fact al¬ 
ready hold that information. 

In this paper we present an analysis of emission 

for a sample of 16 embedded protostars from previously pub¬ 
lished/archival Submillimeter Array observations with the aim 
of revealing their CO emission distributions. The paper is laid 
out as follows: ^describes the sample of sources and presents 
an overview of the observed images and spectra, ^presents an 
analysis of the observed maps in contexts of detailed dust and 
line radiative transfer models - with a specific case study of the 
Class 0 protostar IRAS 03282-1-3035 setting up a more general 
analysis for the sample of 16 sources. The key result, that half of 
these protostars in fact show more extended CO emission than 
what one should expect based on their current luminosities, is 
discussed in along with the statistical implications for the 
frequency of bursts and the properties of the CO ices, ^sum¬ 
marizes the main findings of the paper. 


2. Data 


As basis for this analysis we constructed a sample of embedded 
protostars observed using the Submillimeter Array (SMA; |Ho| 
et al.|20()4] l. One of the preferred spectral setups for such obser¬ 
vations at 230 GHz (1.3 mm) is a setting covering the 7 = 2-1 
transitions of '^CO, ^^CO and C'^O. For the purpose of this pa¬ 
per we focus on C'*0, which typically shows centrally concen¬ 
trated emission toward the location of the protostars. We utilized 
observations in the SMA’s compact-North and compact config¬ 
urations that typically result in an angular resolution of 2-3" 
at these frequencies. Due to the interferometer’s lack of short- 
spacings, the observations are typically not sensitive to emission 
extended over scales larger than about 15". 

The sample was put together utilizing the PROSAC study of 
Class 0 protostars pprgensen et al.|2007| l as well as the sources 
from the SMA archive that were also used as a basis for the anal¬ 
ysis of binarity by [Chen et ffil] ( |2013| l. This sample comprises 
predominantly Class 0 objects - with a few young Class I ob¬ 
jects (i.e., sources with Tboi ^ 150 K - or envelope masses of a 
fewxO. 1 Mq to a few Mq). We exclude sources with low signal- 
to-noise detections either due to their intrinsic C'^O 2-1 line 
strengths or poorer data quality. 

A few sources for which observations have not previously 
been presented were added. Those include IRAS 03256-3055 
observed on 2011/09/24 as part of a larger program to survey 
the population of young stars in the NGC 1333 cluster (PI: 
J. Williams) as well as observations of two intermediate mass 
protostars in Orion, MMS6 and MMS9 ( Johnstone et al.||2003[ 
jjprgensen et al.|[2006b| ), from observations on 2007/02/07 (PI: 
J. Jprgensen). The full sample of sources is summarized in Ta¬ 
ble [U 

Both the archival data and the new observations were re¬ 
reduced using the MIR package 20121 and imaged and 
cleaned using Miriad (Sault et al. |1995|l. In the reduction we 


followed the standard procedures, calibrating the complex gains 
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Table 1. Sample of sources 


Source 

anooo 

bj2000 

J ^ 

^cur 

Distance SMA reference* 

Luminosity reference 


(hh mm ss.s) 

(dd mm ss) 

(Le) 

(pc) 



L1448I2 

03 25 22.4 

+30 45 13 

1.7 

235 

Chen et al. 

(20i: 

5i 

L1448I3 

03 25 36.3 

+30 45 15 

8.8 

235 

Chen et al. 

(20L 

L 

L1448C 

03 25 38.9 

+30 44 05 

8.4 

235 

Jprgensen e 

al. ( 

2007: 

IRAS03256 

03 28 43.4 

+31 17 37 

1.7 

235 

11/09/24 (J 

Wili 

tarns) 

IRAS2A 

03 28 55.6 

+31 14 37 

20 

235 

Jprgensen et al. ( 

20071 

SVS13A 

03 29 03.7 

+31 16 07 

34 

235 

Chen et al. 

(^201: 


IRAS4A 

03 29 10.5 

+31 13 32 

9.9 

235 

Jprgensen et al. ( 

2007 1 

IRAS4B 

03 29 12.0 

+31 13 08 

4.4 

235 

Jprgensen et al. ( 

2007 1 

IRAS03282 

03 31 20.9 

+30 45 30 

1.1 

235 

Chen et al. 

rS- 

li - 

TMRl 

04 39 13.9 

+25 53 21 

3.8 

140 

10/11/29 (F 

uw 

Yen) 

L1527 

04 39 53.9 

+26 03 10 

1.9 

140 

Jprgensen et al. ( 

20071 

MMS6 

MMS9 

05 35 26.0 

05 35 23.4 

-05 05 43 
-05 01 31 

27 

6.3 

420 

420 

U//UZ/Z4 (J. J0rgensen) 
07/02/24 (J. J0rgensen) 

IRAS 15398 

15 43 02.2 

-34 09 07 

1.6 

155 

Chen et al. 

(201311 

CrA32 

19 02 58.7 

-37 07 36 

1.3 

130 

Feterson et al. 

(I201L 

B335 

19 37 00.9 

+07 34 10 

1.3 

100 

Jprgensen et a 

CpD^i 


Evans et al. 

2009 


Evans et al. 

2009 


Green et al. 

(2013 


Evans et al. 

2009 


Jprgensen ef 

al. 

121 

107 

Evans et al. 

20 

19 


Karska et al 

( 

2(jn 

1 

Karska et al 

( 

2013 

1 

Evans et al. 

'z 

C 

C 

Karska et al 

( 

20T3 

1 

Karska et al 

( 

2013 

1 


(see below) 


(20131 


Karska et al _ 

Feterson et all (|2 UTT1 i 

Launhardt et ak ^2(JT3 1 


Notes. Bolometric luminosities compiled from recent Herschel and Spitzer surveys with references in the final column. For the two sources in 
Orion, MMS6 and MMS9, some differences are found in literature. We adopt the value of 27 Lq for MM S6 from|Manoj et al.| ( [2013| l and derive 
the internal luminosity of MMS9 by scaling by the ratio of the 70 um fluxes of the tw o sources reported by |Billot et al.| ( |2012^ (see also discussion 
of the relation between 70 ^m flux and internal luminosity in Dunham et al.|20d8 '. Published paper first describing data - or, if no previous 
paper, date of observations and PI. 


through observations of nearby quasars (typically two per track), 
passband calibration through observations of strong quasars be¬ 
fore and/or after the observations and flux calibrations utilizing 
planet observations. Continuum sensitivities (RMS) in 2 GHz 
bandwidth for these observations were 0.5-1 mJy. Typically the 
C'^O 2-1 observations had been performed by allocating 128- 
1024 channels to one of the (then) twenty-four 104 MHz chunks 
of the SMA correlator placed to cover the transition. This 
resulted in spectral resolutions of 0.14-1.1 km s ' with a typical 
line RMS noise of 30-60 mJy beam * averaged over 1 km s * 
channels. 


and line emission in Fig. for those sources it is seen that the 
line emission is generally better aligned with location of the cen¬ 
tral heating source as traced by shorter wavelength (e.g. Spitzer 
data) and the SMA continuum emission is more strongly affected 
by the dust density and temperature distribution on the scales of 
the envelope-disk transition. 

In the following we focus on reproducing the extent of the 
integrated emission from the C'^O 2-1 transition - but further 
restrict the analysis beyond what is shown in Figs.[T]and[^and do 
not include the most extended (low surface brightness) emission 
as justified below. 


2.1. Overview of data 

Figs. 0 andj^present overviews of the observed C'^O 2-1 spectra 
and maps. The spectra were extracted in the beam toward the 
peak of the continuum emission for each source. The maps were 
made by integrating the emission over +1.5 kms * around the 
systemic velocity for each line. 

The spectra show fairly regular profiles with little evidence 
for wings or other asymmetries as often is the case for the more 
common isotopologues - as well as in larger scale single-dish 
maps. For the maps the emission is also relatively centrally con¬ 
densed with low surface brightness extended emission in a few 
cases - but interestingly not for the most prominent outflow 
sources. Fig. S also shows that for the sources with some struc¬ 
ture in the 2-1 emission there is no clear correlation with 
the outflow axis probed by *^CO 2-1. These observations likely 
reflect that the optically thin integrated C'^O 2-1 emission is rel¬ 
atively little affected by the outflows present on larger scales but 
rather associated with the denser parts of the envelopes on scales 
< 1000 AU - in agreement with what has been inferred from 
previous high angular resolution millimeter wavelength studies 
1 12006 1. In a few cases (e.g., IRAS03256 


(e.g., Arce & Sargent 


and IRAS 15398) sma: 


3. Analysis 

The main aspect of the analysis of the C'^O 2-1 emission is to 
test whether its extent can be reproduced through simple radia¬ 
tive transfer models of the envelope structure and protostellar 
luminosities. We first present an in-depth analysis of a single 
source as a test case ( ^3.1| l before considering the entire sample 
of sources in a general manner (|3.2|l. 


1 offsets are seen between the continuum 


3. 1. Test case: IRAS 03282+3035 

As a test case we focus on the IRAS 03282+3035 Class 0 pro¬ 
tostar (IRAS03282 in the following): IRAS03282 is a deeply 
embedded protostar in the Perseus molecular cloud driving a 
northwest/southeast oriented outflow seen in high resolution mil¬ 
limeter interferometric maps (e.g., Arce & Sargent 2006| l and 
in mid-infrared images from the Spitzer Space Telescope (e.g., 
Jprgensen et al.|2006a i. Its compact C'^O 2-1 emission is par¬ 
ticularly bright and regular, making it a good test case for the 
method. 

To model the C'^O 2-1 profile we construct a model for the 
envelope based on its submillimeter continuum and line emis¬ 
sion - similar to the approach adopted in, e.g., |j0rgensen et al.| 
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Fig. 1. C**0 2-1 spectra toward the peak of its 
integrated emission for each of the protostars in 
the sample. 
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Fig. 2. Integrated C**0 2-1 emission for each 
of the protostars shown as contours on top of 
continuum maps of each source. The contours 
are shown at steps corresponding to 20%, 30%, 
40% ... of the peak emission in the integrated 
line maps. The dashed lines indicate the out¬ 
flow axis from SMA *^CO 2-1 maps for the 
sources where this line was observed simulta¬ 
neously with C'*0 2-1. 
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( 2002|[2009| ). We adopt a 1-dimensional power-law density pro¬ 
file for the envelope, p oc r^P, and calculate its temperature 
profile self-consistently based on the luminosity of the source. 
For the calculations we utilize the radiation transfer code Tran- 
spfier^that solves the ID dust radiative transfer problem for 
absorption and re-emission following the methods described in 
|Dullemond et al.| ( |2002[ ). The normalization of the envelope den¬ 
sity profile is fixed by comparison to single-dish submillime¬ 


ter continuum observations (for further discussion see, e.g., J0r- 


Igensen et al.||2009| l. In the following we test models that fix p 
at 1.5 (the expectation for a free-falling envelope) and 2.0 (cor¬ 
responding to a static isothermal sphere). From radiative trans¬ 
fer models of submillimeter continuum brightness profiles from 
JCMT/SCUBA observations, p is typically found to vary be¬ 
tween these values with uncertainties of +0.2 in the measure¬ 
ments for individual sources (e.g., |j0rgensen et al.|2002[ [Shirley | 
|et al.|200^|Kristensen et al.|2012|l. 

The resulting density and temperature profiles from the dust 
radiative transfer calculations are then used for line radiative 


transfer calculations using the Ratran code ( Hogerheijde & van] 
der Tak|2000|l. Adding the velocity field and an abundance pro¬ 


file for a given molecule, Ratran solves the full non-LTE radia¬ 
tion problem to estimate its level populations as function of lo¬ 
cation and subsequently ray-traces the solution to calculate syn¬ 
thetic images that can be compared directly to the observations. 
For the CO abundance profile we adopt a simple step function 
(e.g. J0rgensen et al.|2005] l. We multiply the synthesized images 
with the interferometric primary beam and Fourier transform the 
images for the direct comparison to the observations in the (m, v)- 
plane. Fig. [^explores the influence of different models on the CO 
emission profiles. In the panels of this figure the data and mod¬ 
els are represented by the amplitude of the integrated emission 
as function of projected baseline length. 


With these choices the model has a number of free parame¬ 
ters; the envelope inner and outer radius, its density normaliza¬ 
tion and the CO abundance profile as function of radius. By re¬ 
stricting the comparison to specific baselines the inner and outer 
envelope radii are not important for the results and we adopt 
25 AU and 8000 AU as in J0rgensen et al.|p009| l. The normal¬ 
ization of the density profile and the absolute abundances are 
naturally degenerate - in the sense that a denser/more massive 
envelope can fit the data with a lower abundance. The shape of 
the abundance profile on these scales is critical, however. 


Fig. [^compares the observed C'^O visibility amplitudes to 
the predictions from a range of models summarized in Table 
Fig. illustrates the simplest models where the C**0 abun¬ 
dance is kept constant throughout the envelope at 2 x 10“^, cor¬ 
responding to a canonical CO abundance of 10 - as well as 

models where the abundance is increased and decreased by a 
factor 8. Two important conclusions can immediately be drawn 
from this plot: if the C'^O abundance is constant at the canonical 
level the emission of the 2-1 transition becomes optically thick 
throughout the envelope - so that further increases in abundance 
do not change the brightness profiles significantly. For low abun¬ 
dances, the model underproduces the emission on large scales 
(short baselines) as well as in single-dish observations. For all 
three models little emission is seen on long baselines compared 
to the observations that show significant emission out to base¬ 
lines of approximately 35^0 kA. This indicates that the abun¬ 
dance must change on small scales - rather than staying con¬ 


* http://WWW.ita.uni-heidelberg.de/ 

~dullemond/software/transphere/index.shtml 


stant throughout the envelope - making the emission maps ap¬ 
pear more centrally condensed. 


Models with an abundance enhancement of CO by up to two 
orders of magnitude (Fig. |^) at the radius where the temper¬ 
ature increases to either 20 or 30 K (corresponding to the CO 
sublimation temperature of pure or slightly mixed CO ice) do a 
slightly better job in the sense of producing some emission at the 
smaller scales/longer baselines - but are still not able to capture 
the full emission profiles. Also, for a canonical CO abundance 
on small scales the model severely underproduces the emission 
on larger scales, so the change has to occur on larger scales than 
what would be predicted given the current luminosities. The fig¬ 
ure also includes a model with an abundance decreasing through 
the protostellar envelope from the outer edge towards the cen¬ 
tre until the point where the the temperature reaches 20 K and 
then increases abruptly back to the canonical value. Such models 
could capture, e.g., the gradual freeze-out with increasing den¬ 
sities ( J0rgensen et al.|2005[ Fuente et al.|2012| inferred on the 
basis of single-dish and (lower resolution) interferometric data. 
In addition, this model does not affect the emission on the longer 
baselines. Thus, the key point from this analysis is that the shape 
of the visibility amplitude curve as function of projected base¬ 
lines is not affected by the exact abundance structure on larger 
scales. The task is therefore to determine what luminosity is re¬ 
quired to shift the CO sublimation temperature out to the radius 
where the enhancement is seen in the CO emission maps. 


Fig. 15 demonstrates three models that fit the data well: these 
represent models where the CO sublimates with an abundance 
enhancement within the radius where the temperature increases 
above 20 K in a p cx r * ^ envelope around a 5 L© protostar as 
well as two models where CO sublimates at 30 K and the density 
profile has either p = 1.5 with a luminosity of 30 L^or p - 2.0 
and a luminosity of 50 Lq. These models are indistinguishable 
at baselines longer than 15 kd. Increasing or decreasing the lu¬ 
minosity by more than about 20% makes the fits significantly 
worse. These models illustrate some of the degeneracy between 
the best fit parameters and the limitations of this method: nat¬ 
urally there is a degeneracy between the temperature at which 
CO sublimates and the required luminosity enhancement - al¬ 
though this can be somewhat alleviated by considering the sam¬ 
ple in its entirety. Likewise the exact luminosity enhancement is 
dependent on the steepness of the density profile: a steeper den¬ 
sity profile causes the C'^O brightness profile to become steeper 
and thereby requires a larger enhancement in luminosity to shift 
the sublimation radius outwards. In the following we continue to 
consider the case where p - 1.5- but note that the required lu¬ 
minosity enhancement may be up to a factor 2 larger for steeper 
density profiles (i.e., in that context the inferred enhancement is 
a lower limit). It is clear that as long as only the longer baselines 
are considered the structure and/or abundance on larger scales 
(shorter baselines) are unimportant for the fits. This point is fur¬ 
ther illustrated in Fig. 15 where the strength of the external inter¬ 
stellar radiation field is varied rather than the internal luminosity 
of the source: such models may vary the profiles on the large 
scales through the sublimation of the CO and changing the exci¬ 
tation temperature - but do not affect the shape of the emission 
profile in the inner 10" (baselines longer than 15 kd) - corre¬ 
sponding to radii <1000 AU for the studied sources. 


In summary, for the test case of IRAS03282 it is not possi¬ 
ble to reproduce the observed C'*0 brightness profiles unless its 
abundance profile reflects sublimation at a temperature set by a 
higher luminosity of the source than its current ^ 1 Lq. 
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Fig. 3. Comparison between different parameterisations of the envelope physical and chemical structure (for details of the parameters for each 
model, see Table comparing the observed and modeled visibility amplitude as function of projected baseline length for IRAS03282. In each 
panel the observed integrated C**0 emission is seen as datapoints with models as lines. The grey histogram indicates the expected amplitude for 
zero signal. In panel a) models with constant abundances and a luminosity corresponding to the current source luminosity (1.1 Lq). In panel b) 
models with a changing abundance at the radius where the temperature increases above 20 K. In panel c) the source luminosity is increased and 
CO consequently sublimates at larger radii providing the variation in intensity as function of projected baselines required to fit the data. Finally 
panel d) tests the importance of varying the strength of the external interstellar radiation field - which again is found to only affect the outer parts 
of the envelopes. 


3.2. Sample approach 


Inspired by the discussion above we define the approach for the 
entire sample of sources. For each source we derive the extent of 
the integrated C'^O 2-1 emission by fitting 2D Gaussian profiles 
in the («, u)-plane to limit systematics due to the deconvolution. 
We restrict the fit to baselines longer than 15 kd (about 20 m; 
or scales corresponding to about 10")- As shown above ( |3.1| l 
this selection is ideal for measuring the extent of the regions of 
the envelopes where CO sublimates - and furthermore filtering 
the remainder of the low surface brightness extended emission 
associated with the ambient cloud. Modeling the emission on 
shorter baselines can be used together with single-dish observa¬ 
tions to probe the amount of CO in the larger scale envelope and, 
for example, the rate at which CO freezes-out or the importance 
of photo-desorption (e.g., Jprgensen et al. [2005 [ Fuente et al. 
|2012| l. However, as we here are mainly interested in the gradi¬ 
ents in C'^O emission associated with CO sublimating at higher 
temperatures and aim to treat all sources consistently, we focus 
only on the emission on the longer baselines. Table lists the 
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deconvolved extent of the C'^O 2-1 emission (FWHM) for each 
source and Fig. plots the extents (as the radius = FWHM/2) 
versus the current luminosities of the sources. Already from this, 
it can be seen that there is not a correlation between the decon¬ 
volved extents and the current luminosities of the protostars - 
what one otherwise would expect if the heating by the central 
protostar dominates the temperature profile and thus determines 
the size of the region where CO is in the gas-phase. 

Fig.j^also shows the predictions for the C^*0 emission from 
a set of radiative transfer models. We do not construct models for 
each individual source - but rather make predictions of the C^^O 
emission extent from a generic envelope model with p oc 
and ranges of luminosities. For each of these models we follow 
the same procedure as described in eu performing the dust and 
line radiative transfer calculations and predicting the C^^O 2-1 
emission profile assuming that the CO comes off dust grains with 
a two orders of magnitude change in abundance for a sublimation 
temperature of 30 K (see the discussion about a possible lower 
sublimation temperature below). The models are processed in 
the same way described above, multiplied by the interferometric 
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Table 2. Summary of models from Fig.[^ 


Panel 

Model 

Density 



Abundance 

Luminosity 

a) 

Al 

p OC 



Xq (constant) 

1.1 Lo 


A2 

p OC 



Ao/8 (constant) 

1.1 Lo 


A3 

p OC 



8 Aq (constant) 

1.1 To 

b) 

B1 

P OC 

Xq for r 

< r20K 

Ao/lOOforr > rioK 

1.1 To 


B2 

p OC 

Xq for r 

< r20K 

Ao/10 forr > rioK 

1.1 To 


B3 

p OC 

10 Aq for r 

< r20K 

Ao/10 forr > rioK 

1.1 To 


B4 

p OC 

Xq for r < 

r 20 K - and gradual decrease from Xq at rout to Aq/IOO at r 20 K 

1.1 To 

c) 

Cl 

p OC 

Xq for r 

< ^20K 

Ao/10 forr > rioK 

5 To 


C2 

p OC 

Xq for r 

< ^'SOK 

Ao/10 forr > rsoK 

30 To 


C3 

-2 0 

peer 

Xq for r 

< ^'SOK 

Ao/10 forr > rsoK 

50 To 

d) 

D1 

p OC r~^'^ 

Xq for r 

< r2QK 

Ao/lOOforr > rioK 

1.1 To 


D2 

P OC 

Xq for r 

< r20K 

Ao/lOOforr > rioK 

1.1 Lo ;ISRFx5“ 


D3 

p OC 

Xq for r 

< r20K 

Ao/100 forr > r 20 K 

1.1 Lo ;ISRFx50'' 

Jotes. Aq 

is a C'*0 abundance of 2 x 10 ’ with respect to H 2 , corresponding to a canonical CO abundance of 10 

Interstellar radiation fiel 


primary beam field of view, Fourier-transformed and the result¬ 
ing extent measured by fitting a 2D Gaussian-profile in the {u, u)- 
plane. The model predictions are compared to the observational 
data in Fig.|^with Table [^providing the current luminosities, the 
predicted sizes of the CO emission profiles given these luminosi¬ 
ties, the actual measured sizes and the corresponding enhance¬ 
ments in luminosities. The uncertainty of the inferred value for 
the luminosity enhancement from the model hts is about 20% 
( ^3.1[ ). Due to systematic uncertainties (e.g., due to the model 
assumptions such as the envelope mass in the general approach), 
variations should only be considered signihcant on a 50% or 
higher level - but generally considered lower limits due to the 
method being insensitive to the most intense bursts due to the 
interferometer resolving out more extended emission and the 
degeneracy with the steepness of the envelope density prohles 
(®. 


A few things are immediately apparent from Fig. the ob¬ 
served extents of the CO emission do not show a clear corre¬ 
lation with the current luminosities with “typical” variations in 
the deconvolved sizes for any given luminosity of about a factor 
of three. The exception to this statement is that no sources with 
high luminosities have very compact C'*0 2-1 emission. If the 
extent of the C^^O emission is determined primarily by the sub¬ 
limation due to heating by the central protostar this is unsurpris¬ 
ing, however: as the sublimation is expected to be instantaneous 
compared to the evolutionary time-scales of the protostars - it is 
unlikely that one encounters such sources. It is indeed seen that 
the sources with the most compact emission for each luminosity 
nicely track the predictions from the radiative transfer models. 
In contrast, for about half the sources the CO emission is seen to 
be extended as one would expect if the luminosity had been from 
a factor about hve up to two orders of magnitude higher than it 
is now. Under the same interpretation this is an indication that 
their current luminosities are not what determines the observed 
extent of the CO emission. 



Lum [Lq] 

Fig. 4. Radii of C**0 emission at the half-maximum intensity for the 16 
Class 0 protostellar sources studied in this paper (symbols) as function 
of luminosity compared to the predictions from radiative transfer mod¬ 
els: the solid line indicates a 1 Mq envelope with sublimation of CO and 
resulting increase in abundance by two orders of magnitude within the 
radius where the temperature increases above 30 K. The shaded area 
corresponds to envelope masses ranging from 0.3-3 Mg. The dotted line 
represents models where the sublimation takes place at a temperature 
of 20 K instead. Finally, the dashed lines correspond to sources where 
the luminosity has been increased by one and two orders of magnitude 
and the sublimation takes place for a temperature of 30 K. The typical 
uncertainties on the measured extents are 10-20% of the plotted values; 
the typical uncertainties in the predicted extents are up to 50% due to 
the uncertainties in the underlying model, luminosity and similar. 

4. Discussion 

4.1. Envelope mass and CO sublimation temperature 

Two parameters that may affect the exact value of the inferred 
increase in luminosity on the basis of the CO extent are the as¬ 
sumed envelope mass and CO sublimation temperature in the 
generic models. 

A simple test reveals that the former is of little importance 
for the solar-type protostars considered here: Fig. shows the 
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Table 3. Measured and predicted extents of the deconvolved C'*0 emission - as well as the current luminosities and luminosity enhancements 
required to reproduce the extended C'*0 emission. 


Source 

^cur (^ 0 ) 

Predicted size (")“ 

Measured size (")* 

L'CO / ^cur 

L1448I2 

1.7 

0.9 

5.3 

25 

L1448I3 

8.8 

2.2 

4.9 

4.2 

L1448C 

8.4 

2.2 

2.0 

0.8 

IRAS03256 

1.7 

0.9 

3.4 

12 

IRAS2A 

20 

3.3 

3.0 

0.8 

SVS13A 

34 

4.6 

6.7 

2.0 

IRAS4A 

9.9 

2.4 

5.5 

4.6 

IRAS4B 

4.4 

1.6 

1.9 

1.4 

IRAS03282 

1.1 

0.8 

4.0 

25 

TMRl 

3.8 

2.4 

5.7 

5.4 

LI 527 

1.9 

1.6 

2.6 

2.3 

MMS6 

27 

2.0 

4.1 

3.7 

MMS9 

6.3 

1.1 

3.6 

9.6 

IRAS 15398 

1.6 

1.3 

3.5 

5.9 

CrA32 

1.3 

1.5 

7.9 

24 

B335 

1.3 

1.9 

3.2 

2.7 


Notes. Predicted extent of the C**0 2-1 emission given the current luminosity of the source (Lcur)- ® Measured deconvolved extent of the C**0 
2-1 emission (FWHM). The uncertainties are dictated by the signal-to-noise and the assumption of the emission being distributed as Gaussians: 
typical statistical uncertainties are 0.1-0.5" - or 10-20% of the measured sizes. Ratio of luminosity needed to reproduce C'*0 2-1 extent to 
current luminosity - i.e., inferred enhancement in luminosity. See the text for discussion about the assocated uncertainties. 


spread in inferred CO sublimation radii for envelope masses 
varying by an order of magnitude from 0.3 to 3 Mq: the vari¬ 
ation in the predicted CO sublimation radii is about 10%. This 
reflects that the envelopes are largely optically thin to their own 
radiation in which case the temperature profile can be described 
by a power-law temperature profile with radius with the over¬ 
all scaling only determined by the luminosity of the source (e.g., 


Scoville & Kwan|1976[ Doty & Leung 1994[ Chandler & Richer 


20001 1 . 


The adopted CO sublimation temperature is much more im¬ 
portant as also illustrated in Fig.|^ Shifting the sublimation tem¬ 
perature from 30 K to 20 K in the models shifts the radius where 
CO evaporates out to align with the predictions for where it 
would be for a model with a sublimation temperature of 30 K 
but luminosities an order of magnitude higher. 

This discussion is highly relevant: traditionally many stud¬ 
ies of modeling sublimation of CO in star-forming regions adopt 
a binding energy of 855-960 K from laboratory experiments of 
pure CO (e.g., Sandford & Allamandola|[l993[ Bisschop et al. 
[2006 | l. However, laboratory experiments also show that mix¬ 
ing the CO with H 2 O ice can increase the CO binding en¬ 
ergy to ai 1200-1300 K (e.g., |Collings et aT]|2003[ [Noble et al. 
|2012^ . The thermal sublimation rate at a constant dust temper¬ 
ature, Td, depends exponentially on the binding energy, Et,, as 
^ oc exp{Et,ITd) (e.g., Eq. 4 of [Rodgers & Charnley|2003] l. For a 
binding energy of 960 K it takes less than 1 year for CO to sub¬ 
limate for a dust temperture above 21 K. For a binding energy 
of 1300 K the sublimation rates slows down so that it takes less 
than 1 year for CO to sublimate only at temperatures above 28 K. 
Observationally, there is some evidence for a higher sublimation 
temperature based on modeling of multi-transition CO observa¬ 
tions as well as previous high angular resolution observations 


( j0rgensen[2OO4t Jprgensen et al.|2005||Yildiz et al.|2()T3 1 . 

It is worth noting that half of the sources in fact do show 
compact CO emission that is well-fit by the models with a subli¬ 
mation temperature of 30 K adopting their current source lumi¬ 
nosities. That directly argues against a lower sublimation tem¬ 
perature as it would result in more extended CO emission than 
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observed. In contrast, a higher sublimation temperature cannot 
be ruled-out based on the observations here - but it would im¬ 
ply a larger fraction of the sources have experienced increases in 
luminosities in their recent histories. 


This result is not necessarily a contradiction with some high 
resolution imaging studies of more evolved disks around T Tauri 
stars that may favor lower sublimation temperatures through in¬ 
direct imaging (e.g., IQi et al.||2013| [Mathews et al.||20T3] l - al¬ 
though it might be interesting to revisit those studies with mod¬ 
els for which a higher binding energy/sublimation temperature 
would be adopted. Whereas the CO and H 2 O ices are likely 
formed simultaneously in the protostellar environments where 
the temperatures generally are low - the ices in the disks around 
T Tauri stars may have undergone some heating and sublima¬ 
tion before reformation when entering the disks and the material 
cooling down. As the binding energy of H 2 O is much higher than 
that of CO, only a small fraction of the water may not sublimate 
before entering the disk (e.g., Visser et al.[[2009[ [Cleeves et al 


2014|l. This would naturally lead to the formation of a clearer 


onion-shell structure of the ices in the disks with an inner layer 
of H 2 O ice and an outer layer of more purified CO ice with a re¬ 
sulting lower binding energy. It is possible that the CO sublima¬ 
tion temperature varies from source-to-source, e.g., depending 
on the water content in their ambient environment. This would 
naturally be an important result in itself - although, no other 
direct observational evidence exists for such variations. In that 
case, one could only argue that three out of the sixteen sources 
would show evidence for enhanced luminosities in their recent 
histories. 


In any case, this discussion emphasizes the need for high 
quality laboratory measurements of ices relevant for astrophysi- 
cal environments - and a careful consideration of the exact con¬ 
ditions and chemical histories of the studied regions. 
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pre-burst 

t > fdep since last burst 


burst 



Fig. 5. Cartoon illustrating the interpretation 
of the extended CO emission signatures due 
to variations in source luminosities. The dark 
blue and yellow regions indicate where CO is 
frozen-out and in the gas-phase in three char¬ 
acteristic situations around a protostellar burst 
in luminosity. The red dashed line is the radius 
where the temperature in the envelope reaches 
the sublimation temperature for CO. The three 
sketches indicate the situations (1) just before a 
burst where the CO sublimates out to a radius 
defined by the current luminosity of the proto¬ 
star, (2) during a hurst where the luminosity of 
the protostar increases and the CO sublimation 
radius moves outwards, and (3) immediately af¬ 
ter a burst where CO remains in the gas-phase 
in a large region - before it freezes out again. 


4.2. Ice sublimation as a tracer of accretion bursts 

A fairly straightforward interpretation of the variations in the ex¬ 
tent of the central C'^O emission relative to the predictions given 
the current source luminosities is that the protostars in fact un¬ 
dergo bursts of accretion and corresponding luminosities. Fig.[^ 
illustrates this progression through three characteristic phases: 
in a period of time before the burst, or equivalently a long time 
since the previous burst, the chemistry will set itself in a rela¬ 
tively steady state with CO frozen-out outside the radius where 
the temperature reaches the sublimation temperature and in the 
gas-phase within this (Stage 1 in Fig. [^. During the burst this 
radius will be shifted outwards and due to the short time-scale 
for sublimation of CO the freeze-out/sublimation boundary will 
shift outwards correspondingly (Stage 2). For typical binding en¬ 
ergies for CO (see ^4.1| l the time-scale for CO to sublimate at 
temperatures of 30 K is of order minutes to days - i.e., much 
shorter than the typical dynamical time-scales related to the in¬ 
fall of material from the envelope. Finally, after the burst has 
finished and the source reset to a low luminosity, CO will begin 
to freeze-out again. However, due to the slower freeze-out it will 
still appear to be extended for a longer period of time (Stage 3). 

It is possible to be more quantitative about the time-scales. 
For a characteristic envelope with a mass of 0.3-0.5 Mq, a size 
of 10,000 AU (radius) and a density profile dropping as p cx 
r ' the time-scale for CO to freeze-out at the den sity at 300- 


500 AU is of order Uep ~ 10^ years (e.g., Eq. 3 of Rodgers & 


Charnley||2003|l. Thus, if a given source undergoes an increase 


in luminosity, CO would stay in the gas-phase - and its emission 
appear extended - for this period of time after the luminosity has 
decreased again. 

With additional assumptions one can even turn this into a 
statement about the frequencies of such bursts implied by the 
current data. One can for example assume that the duration of 
the burst in Stage 2 ( 12 ) from Fig. is much shorter than the 
total duration of Stage 1 through 3 (fi + t 2 + h): based on the 
observed luminosity distributions of young stellar objects from 
the Spitzer Space Telescope c2d legacy program and simple cal¬ 
culations of the accretion rates, [Evans et al.| ( |2009| ) argued that 
half the final mass of a given young star is accreted during 7% of 
the duration of the embedded stages. Similarly from a compari¬ 
son between numerical hydrodynamical simulations of collaps¬ 


ing cores and the distributions populations of young stars, Dun 
|ham & Vorobyov| ( |2012| l argued that modeled protostars on av¬ 
erage spend ^1% of their lifetimes in modes of accretion bursts. 
The analysis above shows that half of the sources have CO ex¬ 
tents corresponding to luminosities about a factor five or more 


above the current. In context of the scenario outlined above this 
would imply that fi ^ t^. As the duration of per definition is the 
depletion time, Uep, this suggest that a given source undergoes a 
burst every fi H- f 3 ~ 2f(iep ~ 20,000 years - or of order 5 bursts 
during the 10^ year duration of the deeply embedded Class 0 
stage. If one on the other hand assumes that variations in the 
CO sublimation temperatures between the different sources is 
the dominant reason for the variations in the CO extent and thus 
that only three out of the sixteen sources show significant C^^O 
2-1 extents, this would imply a burst every 50,000 years. These 
time-scales are in good agreement with estimates by |Scholz et al.| 
P013| l who compared Spitzer and WISE photometric data for 
young stars in nearby star forming regions and identified 1^ 
bursts in a sample of 4000 young stars. Translated into frequen¬ 
cies these numbsers would correspond to bursts in intervals of 
about 5000 and 50,000 yr. 

It is worth re-emphasizing the reasons why the method out¬ 
lined in this paper is useful: first and foremost the C'^O emission 
is not significantly affected by, e.g., outflows or the ambient en¬ 
vironment. This is observationally shown given the narrow lines 
and relatively centrally condensed distributions of the line emis¬ 
sion. This is mainly a result of the C'^O being largely optically 
thin and thus predominantly sensitive to the high column den¬ 
sities of material associated with the protostars themselves. An 
analysis of the chemistry of common molecular tracers during 


and after the accretion burst is presented by Visser et al. ( 2015| l. 
Also, any more smoothly distributed material on larger scales in 
the clouds is filtered by the interferometer. Secondly, the time- 
scales for sublimation and freeze-out of the CO molecule are 
apparently well-matched for its use as a tracer. Observationally 
the typical envelope masses (and consequently densities) change 
on time-scales of a few 10^ years from one to a few Mq to 0.1 to 
0.5 Mq for early Class I sources, which is still sufficiently long 
for freeze-out and sublimation to occur. In fact, the observation 
that some sources do show the compact C^^O emission directly 
shows that large variations on short time-scales are unlikely: had 
the frequency of bursts been higher than 1 / fdep the sources would 
never get back to the Stage 1 and cycle between Stages 2 and 3 
- with the result that the CO emission would always appear ex¬ 
tended. 

Of course these estimates are still rather crude due to the 
low number sample statistics. Also, due to the inherent uncer¬ 
tainties in the modeling discussed above, the analysis does not 
constrain the strengths of the bursts. Eor example, due to the low 
sublimation temperatures the CO isotopologues are not sensitive 
to very intense bursts that may cause it to evaporate out to the 
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largest scales probed by the interferometric data (~ 1000 AU 
radius). Observations of other species that sublimate at higher 
temperatures (or species reacting to those) could help refine the 
accretion histories of protostars. Observations of the H 2 O could 
for example be used to directly test whether the binding energy 
of CO varies significantly from source to source: if the extents 
of the H 2 O emission is found to vary in conjunction with CO, 
it would indicate that luminosity variations dominate over varia¬ 
tions in the binding energy. Surveys of larger (unbiased) samples 
of sources coupled to models would greatly improve the statis¬ 
tics also on smaller scales and/or correlate to, e.g., physical prop¬ 
erties of the ambient environment or the physical structure of the 
protostars on small-scales that may drive possible accretion vari¬ 
ations. A more detailed description of the dynamical evolution 
of the protostellar envelopes on the characteristic 10"^ year time- 
scale for freeze-out would also be needed to refine the statistics. 
Finally, although the maps do not show any indications that the 
C'^O 2-1 emission is affected by the presence of outflows, sensi¬ 
tive high angular resolution observations will be needed to quan¬ 
tify what impact outflows have on the quiescent gas and how they 
contribute to shaping the envelopes. 

Turning the argument around: results such as these empha¬ 
size the need to take the accretion histories of embedded proto¬ 
stars into account when considering their chemistry. If similar 
accretion variations can be confirmed to take place during the 
evolution of embedded protostars in general they may have pro¬ 
found implications for the gas-grain chemistry interplay taking 
place during the early stages of protostellar evolution. 

5. Summary 

An analysis of maps of C'^O J = 2-1 emission from the Sub¬ 
millimeter Array (SMA) toward a sample of deeply embedded 
protostars has been presented. The maps trace the distribution 
of C'*0 in the envelope on few hundred to thousand AU scales. 
The maps are compared to detailed continuum and line radiation 
transfer models - in particular, the predictions for the extent of 
the CO emission given the current luminosities of the sources. 
The conclusions are as follows: 

- The C'*0 spectra appear relatively symmetric with little ev¬ 
idence for strong outflow action. The integrated emission 
maps are likewise regular with only faint low surface bright¬ 
ness extended emission. These results indicate that the inte¬ 
grated C'^O emission is predominantly picking up the mate¬ 
rial in the denser parts of the envelopes on < 1000 AU scales. 

- A case study of the IRAS03282 protostar is used to demon¬ 
strate a method in which the deconvolved extent of the C'^O 
emission measured in the (u, t;)-plane can be compared to 
radiation transfer models for the protostellar envelope. The 
key parameter in this comparison is the assumed luminosity 
of the source and consequently the radius at which CO subli¬ 
mates at temperatures of 20-30 K. This dominates over other 
effects such as the envelope physical structure and exact ab¬ 
solute abundances. 

- In the comparison for the entire sample it is found that half of 
the sources show C'*0 line emission that is more extended 
than predicted by models based on the current luminosities 
of the sources if CO sublimates at temperatures of 30 K. For 
those sources the increase has to have been more than about 
a factor five; for four of these the inferred luminosity would 
be more than a factor 10 higher. 

- A CO sublimation temperature lower than 30 K cannot be 
reconciled with the observations of half the sources with 
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compact emission. The sublimation temperature of 30 K in 
these environments may be taken as evidence of a higher 
binding energy for the CO than what is expected for pure 
ices - a possibility if CO is mixed in with CO 2 and/or H 2 O. 

- The half of the sources in the sample that show more ex¬ 
tended C^^O emission are the candidates for those that have 
undergone a recent change in luminosity, e.g., due to a burst 
in accretion rate. This would have had to take place within 
the last ~ 10"^ years, which is the time it takes for CO to 
freeze-out on dust grains at the densities characteristic of 
the 300-500 AU radii of the protostellar envelopes observed 
here. If such bursts of accretion are taking place in general, 
the small number statistics based on this sample would im¬ 
ply that a protostar undergoes of order 5 bursts during its first 
10^ years. If one assumes that the CO sublimation tempera¬ 
ture varies between 20 and 30 K from source-to-source, a 
conservative estimate would be an interval of 50,000 years 
between the bursts. 

This study highlights the importance of studying the phys¬ 
ical and chemical evolution of embedded protostars in unison. 
CO is one of the key species in the gas-phase and grain-surface 
chemistry taking place in the environments of these protostars - 
and understanding the physical mechanisms leading to its freeze- 
out and sublimation therefore of key importance and something 
that needs to be taking into account any models for protostel¬ 
lar chemistry. Conversely, a good understanding of its chemistry 
may be used as a diagnostic tool for effects such as variations in 
the accretion histories of the protostars. An important next step 
would be to extend the current study to a more systematically 
selected sample of sources and differing molecular species, e.g., 
probing different temperature and density regimes, to for exam¬ 
ple break the degeneracies introduced by the limitations in our 
understanding of the details of the chemistry of just one species. 
In the near-future ALMA will make this possible and thereby 
open up the potential for investigating the “dynamic chemistry” 
in star and planet-forming regions. In this context it will provide 
excellent maps for a large number of sources and species and 
provide good statistics also revealing the magnitudes of possible 
accretion bursts. 
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